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(~ l , Context. The Cygnus 0B2 association, the central engine of the Cygnus X star-forming region, is the subject of an extensive 

^ I' INTEGRAL Key Project that will accumulate 6Ms of observations. Analysis of 2Ms of observations by De Becker and co-workers 

provides the most sensitive limit yet obtained on hard X-ray emission from the cluster. 

Aims. We investigate the X-ray emission in the 20-40 keV band expected from the flaring low-mass stellar population in Cygnus 0B2. 
We discuss whether such emission needs to be considered in the interpretation of existing and future X-ray observations of the region, 
and whether such observations might provide insight into the high-energy processes on low-mass pre-main sequence stars. 
Methods. The total hard X-ray flux from low-mass stars is estimated by assuming the observed soft X-ray emission stems from a 
superposition of flares. We further assume the ratio of hard X-ray to soft X-ray emission is described by a scaling found for solar 
flares by Isola and co-workers. 

Results. We estimate the low-mass stellar hard X-ray flux in the 20-40 keV band to lie in the range ~ 2 x 10^' - 6 x 10^' erg .s"' and 
discuss some potential biases that might affect this result. 

Conclusions. Hard X-ray emission could lie at a level not much below the current observed flux upper limits for Cygnus 0B2. If 
this emission could be detected, it would provide insight into the hard X-ray production of large flares on pre-main sequence stars. 
We highlight the penetrating power of hard X-rays from low-mass stellar populations as a possible pointer to our Galaxy's hidden 
star-forming clusters and super-clusters using more sensitive observations from future missions. 
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J>. 1. Introduction mildly relativistic energies (the "thick-targ et model", e.g. lBrownl 

119711; iHudsonll 1 972t iLin & Hudsonll 1 976l) . Tentative direct evi- 
Motivated by the detection of nonthermal radio emission from dence for such nonthermal hard X-rays up to 200 keV has re- 
early-type stars i n the nearby massive Cygnus OB2 associa- cently been found from Swift observations of a large flare o n 
tion_W = 1.7 kpc lMassev & Thompsonlll991h . [be Becker et alj the RS CVn-like interacting binary II Peg lOsten et alj ( l2007h . 
(l2007h have recently used an extensive set of observa- Large stellar flares observed on the young single star AB Dor 
tions obtained by the International Gamma-Ray Laboratory and t he active binaries A lgol and UX Ari by the BeppoSAX mis- 
(INTEGRAL) IBIS instrument to search for hard X-ray emis- sion dBoella et alJI 19971) in the 0.1-100 keV range revealed evi- 
sion from some of the more prominent radio sources. The non- dence of similarly energetic emission, though it was not possible 
thermal radio emission from these stars is thought to arise to ascertain vvhether this was of thernial or n onthermal origin 
from a relativistic population of electrons accelerated by dif- ("M aggio et al.l l2000t iFavata & SchmittI 1 19991: iFranciosini et al.l 
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fusive shock acceleration i n their interacting s tellar winds (e.g. 2001J). 

iPittard & Dougherty! |2006|). |De^Becker et at (2007) note that while any single flare event on a low-mass star is a trifling 
hard X-rays might plausibly be produced by inverse Compton insignificance in comparison with the effusive vigor of its high- 
scattering of photospheric UV photons. mass brethren, the comparative multitude of the former popula- 
One other potential source of hard X-rays in Cygnus OB2 tio" in Cygnus OB2 might render a palpable combined effect. An 
are the myriad flares thought to arise from magnetic recon- estimate of this eff^ect is motivated from diff'erent standpoints, 
nection events that characterize and s ustain the coronae of the Firstly, the search for hard X-rays from early-type stars by 
low-mass stellar population (see, e.g., lAlbacete Colombo et alj iDe Becker et al.l (l2007h resulted in only upper limits; an estimate 
[2007a; Caramazza et al. 2007). The flare paradigm developed of the flux limit at which low-mass stars might be detected could 
from the solar perspective posits that hard X-rays can arise from be important for understanding the origin of any hard X-rays de- 
the impact on dense photospheric gas of electrons accelerated to tected in deeper observations by INTEGRAL or future missions. 
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Secondly, the large concentration of low-mass stars in the 
Cygnus OB2 region might also provide insight into the high- 
energy nature of stellar flares that is diflicult to attain from 
observa tions of nearby s ingle stars that rarely exhibit large 
events. Isola et all (|2007|) have recently shown that a well- 
defined power-law relationship between soft and hard X-rays 
seen in solar flares matches the 20-40 keV flux observed during 
the AB Dor, Algol and UX Ari flares, but t ends to overpredict 
the, pr esumed nonthermal, 60-80 keV flux. [Alexander & Daoul 
(|2007|) find evidence for a saturation in nonthermal X-ray flux 
above 20 keV for solar flares with increasing soft X-ray lumi- 
nosity, and it is conceivable that this is related to the deficit in 
observed nonthermal stellar X-rays. Detection of, or limits on, 
the hard X-ray flux from the low-mass Cygnus OB2 population 
could provide constraints for extrapolating solar flare models to 
much more active pre-main sequence stars. 

Thirdly, the hard X-ray environment of young pro- 
tostellar systems is potentially relevant for the ioniza- 
tion of circumstellar material and protoplane t ary disks (e.g. 
Glassgold et al.' '1997"; 'Aikawa & HerbstI 119991: iFromang et all 
2002; Matsumura & Pudritz 2003). We currently have no infor- 
mation on the intensity of this radiation environment. 

We outline in Sect. |2] below a calculation of the hard X-ray 
flux from the Cygnus OB2 low-mass population under the as- 
sumption that the observed soft X-rays from these stars is dom- 
inated by a superposition of flaring events. We then discuss in 
Sect. |3]the uncertainties in this estimate, the degree to which it 
likely overestimates the true level of hard X-ray emission, and 
the implications of the results. 

2. Estimating the low-mass stellar hard X-ray flux 

There are three main components to our hard X-ray flux es- 
timate: the Cygnus OB2 initial mass function (IMF) and total 
stellar mass: the relationship between stellar mass and observed 
soft X-ray luminosity for the low-mass population: and the rela- 
tionship between the soft and hard X-ray luminosities. In order 
to gain a better perspective on the uncertainties involved in our 
estimate, we use three different approaches that employ these 

ingredients in slightly diff'erent ways. 

We adopt the IMF derived by iKnodlseded (l2000al) based 
on 2MASS photometry. That study noted the power-law in- 
dex found, F - -1.6 + 0.1, is compatible with that of 
iKroupa et al.l (Ll993|) for the solar neighborhood, but is signifi- 
cantly steeper than fou nd from spectroscopy of the central re- 
gion of Cygnus OB2 bv lMassev & Thompson! (1 1 99 lb . We there- 
fore also investigate below the eff'ect of a difl'erent power-law 
slope, and in particula r that found fo r the O rion Nebula Cluster 
(ONC: r = -1.2) by iMuench et all ( 12002.) . We consider hard 
X-rays from stars in the mass range 0.3-2Mo: the upper limit 
corresponds to the approximate boundary between low-mass 
pre-main sequence stars with outer convection zones and in- 
termediate mass stars that should have radiative en velopes: the 
lower mass limit corresponds to the limit of the iKnodlsedeJ 
( l2000ah study — essentially our Umit of current knowledge. For 
stars with mass $ O.3M0, X-ray luminosities of T Tauri stars 
are observed to decline much more rapidly with decreasing 
mass than for higher mass st ars (e.g. iPreibisch et alj 120051: 
lAlbacete Colombo et al.ll2007bl) . rendering their contribution to 
the total X-ray luminosity less significant, regardle ss of whether 
the IMF flattens toward lower masses or not (e.g. IScalollT986t 
iKroupa et al. 1993: Kroupa 2001: Chabrier 2003|). Nevertheless, 
we also investigate the effects of this lower mass cut-off in our 
calculations below. 



lAlbacete Colombo et al.l ( l2007ah found the X-ray luminosi- 
ties of low-mass stars in Cygnus OB 2 to be essentially identical 
to those in the ONC dPreibisch et al.l l2005). We therefore adopt 
the relation between stellar mass, M, and X-ray luminosity, Lx, 
from the ONC study. In order to examine the effects on our cal- 
culation of the significant scatter in observed Lx vs M, we also 
calculate the expected hard X-ray flux using the observed ONC 
stellar sample renormalized to the Cygnus OB2 low-mass con- 
tent (see Sect. 1231). 

The relationship between the soft and hard X-ray luminosi- 
ties was derived using the extrap olation of the relation found 
for solar flares bv llsola et all (120071) . F20-40 = 9.9 x lO^F^", 
where ^20-40 is the 20-40 keV flux density in ph cm"^ keV"' s"' 
and Fc is the Geostationary Operational Environment Satellite 
(GOES) 1 .6-12.4 keV flux in unit of W nT^. In order to convert 
this relation to the stellar case, we assumed that the observed 
soft X-ray luminosities of the Cygnus OB2 low-mass pre-main 
sequence population could be explained by a superposition of 
flares. This assumption is based on evidence that the solar corona 
is largely characterized by a power-law flare dis tribution in total 
energy, E, and frequen c y, A^, dNJdE oc £ " (e.g. lLin et al.lll984l : 
lKrucker&B"en3ll998l: lHudson|[T99l . and on studies of EUV 
and X-ray flares and photon arrival times for active stars that sug- 
gest a similar flare distribution \yith a ~ 2 dAudard et al.ll2000l : 
Kashvap et al. 2002: Giidel et al. 2003). The high-energy tails of 
flares detected in Chandra studies of the ONC and Cygnus OB2 
are also well-approximated by a power-law frequency di stribu- 
tion with a = 2.2 + 0. 2 (ONC, 'Caramazza et al."2007'), a = 
2.1+0.1 (Cygnus OB2. lAlbacele Colombo et al. 2007a). In the 
following we outline our estimates, and then discuss the uncer- 
tainties of this approach. An alternat ive approach t o estimate the 
nonthermal emission is proposed in (lGuedelll2009l) . 

2.1. Hard-soft X-ray luminosity relation for young low-mass 
stars 

The relation from llsola et all (l2007l) relates the flux in the 
1.6 - 12.4 keV band to the flux density in the 20 - 40 keV 
band. In order to obtain a relation between the peak luminos- 
ity of flares in the two bands, we first converted the GOES flux 
from W m"^ into erg s ' luminosity, and the flux density in the 
hard band, given in unit of ph cm"^ keV"' s"', into erg s"' lu- 
minosity, assuming a photon energy of 20 keV and multiplying 
for the band range. We also converted the 1.6 - 12.4 keV lu- 
minosity into luminosity in the 0.5 - 8.0 keV band: the conver- 
sion factor {Kconv - 0.34) was evaluated considering a Raymond 
Smith model wit h kT = 1-35 ke V and logiNn) = 22.25 cm^^ 
( lAlbacete Colombo et alj|2007bl) . We then obtain die following 
relation: 

ZmQ-m = 1-8 • 10''^ ■ (1.2 ■ lO-^'f ■ (i://(o.5-8.0)f erg s'' (1) 

wh ere yS = 1.37 + 0.07 is the index found by llsola et all 
(|2007|) . Note that the factor in the previous relation varies 
strongly due to the uncertainty of p. 

This relation can be converted in a relation between the total 
energy output of flares in the two spectral bands, considering that 
E - X.fi ■ T. Given the power-law distribution of flare energies 
in the 0.5 - 8.0 keV band, we conclude that the 20 - 40 keV 
band flare energies are also distributed as a power-law. The index 
of this power-law, y, is a function of a and /?: y - (a -\- /3 - 
l)/y6, whe re 2.0 < a ^ 2.4 (Caramazza et al. 2007) and 1.30 < 
j3 < 1.44 (llsola et al.ll2007l) . If we impose the condition that y 
also has a value greater than 2, i.e. that the 20 - 40 keV band 
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emission is due entirely to flares, then the acceptable values for 
a and /3 faU in the ranges 1.31 < yS < 1.39 and 2.32 < a < 
2.4. Integrating the distribution of flares between the minimum 
energy and the maximum energy, that we can set as infinite, we 
obtain the following relation between stellar luminosity in the 
20 - 40 keV and 0.5 - 8.0 keV bands: 



L20-40 - A ■ Lo.5-8.0 erg s 



(2) 



where 



/ a-B-1 ^ I ^ ' 

' 0.5-8.0 



Here, we set To.5-8.0 to 10 hours dCaramazza et al. I l2007h . while 
T20-40 is scaled using the ratio (0.16) derived from the me- 
dian duration of flare s in the hard and soft bands calculated by 
IVeronig et al.l (|2002|) from a sample of solar flares. /io.s-s.Omm 
was set to ~ 2 ■ lO-'^ erg that is the median values obtained from a 
statistical analysis of the ONC lightcurves interpreted assuming 
the emission enterely du e to flares with a power law distribution 
dCaramazza et al.|[2007|). The ener gv of flares that contribute to 
the star emission in ICaramazza e t al. (2007) is given in counts. 
We converted the counts to energy, assuming a Raymond-Smith 
emission with the median values of N// and kT of the ONC sam- 
ple (Nh = 1.9 • lO^'cm-^ kT = 1.25 keV). Note that, since /? is 
not too different from 1, the hard x-ray luminosity is relatively 
insensitive to change in -Eo.s-s.O/mn or To.5-8.0 values: a factor of 
10 change in either of these quantities results in only a factor 
of 2 difference in A. Assuming uniform distributions of a and 
/3 in the total intervals given above, we estimate that 90% of 
the resulting values for the proportionality factor A lie in the 
range A=[1.8 ■ 10""^, 3.0 ■ 10"^]. Allowing also for significant un- 
certainties in the other parameters (T20-40/T05-8 o=[0.08, 0.24], 
To.5-8.o=[3, 18] hours, £o.5-8.0m/«=[2- 10^1,2 ■ lO^^] erg), the 90% 
interval for A becomes [1.1 • 10""*, 4.0 ■ 10"^]. We will adopt in 
the following the latter uncertainty range. 



2.2. Analytical Estimate 

Adopting the relation between stellar mass and the X-ray lumi- 
nosity found for the Chandra Orion Ultra-deep Pro ject (COUP) 
stars in the range 0.3 - 2M@ (iPreibisch et al.ll2005h . 



log (Lo.5-8.0) = 30.37 H- 1 .44 ■ log (M) 



(4) 



and combining this with Eqn.|2] we obtain a relation between 
stellar mass and hard X-ray flux. This expression can then be in- 
tegrated over t he stellar mas s distribution for Cygnus 0B2. The 
IMF found by (lKnodlsedej2000a) follows the standard form 



dM 



with F = -1.6, 



(5) 



and assuming for th e present a total mass of M,o, - 5 ■ 10 
(lKnodlsedej2000bh . and maximum and minimum stellar masses 
in the cluster of M^a.^ - 8OM0 and M„„„ = O.3M0, respectively, 
the normalization constant is A: = 1.5 ■ 10"*. 

In order to calculate the total hard X-ray luminosity in the 
cluster, we transform the IMF into a function of 20-40 keV stel- 
lar X-ray luminosity, and integrate between the luminosities cor- 
responding to 0.3 and 2.OM0, obtaining values in the following 
range: 



4o'40 = [l-8-10",6.6-10^2j ^^^^-1 



2.3. Scaling the Orion Hard X-ray luminosity 

We can also make a rough estimate of hard X-ray flux from 
Cygnus OB2 by scaling directly the observed X-ray luminosity 
distribution of the low-mass stars in Orion. This approach has 
the advantage that the observed scatter in the COUP sample is 
intrinsically included. The IMF of the Ori on sample has a shal - 
lower slope tha n that for Cygnus OB2 of lKnodlsede^ (l2000al) . 
with F = -1.2 (iMuench et al.l |2002'), and we investigate the in- 



fluence of this difference in Sect. |2.4| below. 

The total 20-40 keV luminosity for the COUP sample is sim- 
ply the sum over all stars of mass M, in the range O.3-2M0, 

2 

^2o''-40 = ^ ■ Z ^-5-8.0(0 = [8-6 ■ 102^3.1 • 10^"] erg s''. 



where the relation between the 20-40 keV and 0.5-8 keV lu- 
minosities is given by Eqn. (|2]i. Assuming a similar IMF for 
both Orion and Cygnus OB2, the hard X-ray luminosity of the 
latter is then simply given by the product of L^q_^J^ and the 
ratio of cluster masses within the O.3-2M0 mass interval. The 
Cygnus OB2 IMF from Eqn. |5] yields a total mass in this range 
of Mq 3^2 = 3.5 • 1O^M0, while the analogous total mass of the 
COUP sample is M^^^^ = 175M0. Even in this case the total 
Cygnus OB2 hard X-ray luminosity is therefore 



rCyg 
■"20-40 



= [1.7 ■ 10^\6.3 ■ 10"^] ergs 



32^ 



2.4. Simulating Cygnus 0B2 using the Orion sample 

In order to account for the different IMF slopes in Cygnus OB2 
and Orion (F = - 1 .6 c/- 1 .2), we also performed a simple Monte 
Carlo simulation of 10000 stars. The mass range 0.3 - 2 Mq was 
divided into 10 bins, and within each mass bin 0.5 - 8.0 keV 
luminosities were randomly selected from those observed in the 
COUP sample. The number of "stars" drawn in each bin was 
weighted according to the Cygnus OB2 IMF from Eqn. (|5]l, and 
the luminosity for each star was scaled to L20-40 as described 
earlier The luminosity of the simulated sample was then scaled 
so as to have a total mass equal to that of Cygnus OB2 for the 
O.3-2M0 mass range. In this way, we estimate 



r Cyg 
■"20-40 



[1.9- 10^', 6.6 



10^^] ergs-\ 



As a verification of the calculation, we recovered a luminos- 
ity of L^o^f' ^ 9 . 1028,3.2 ■ 10^°] erg i^i— essentially the 
same as that found from direct calculation in Sect. 12.31 — for an 
IMF with slope F = -1.2 scaled to the appropriate COUP total 



3. Discussion 

The range of 20-40 keV fluxes we have evaluated above can 
be compared with the the 3cr 20-60 keV flux upper limit of 
6.1 x\Q-^^ erg crn^ s'^ obtai ned from 2.12 M s of INTEGRAL 
IBIS-ISGRI observations by iDe Becker etaTI (12007 ) for the 
unidentified y-ray source 3EG J20033-i-4188that lies in the close 
vicinity of Cygnus OB2Q. Adopting a distance of 1.7 kpc for 
Cygnus OB2 (e.g. Massex & Thompson 1991.) . this flux corre- 
sponds to a luminosity ~ 2 x 10^^ erg s"'. 



(6) 



' The distance between the position of 3EG J20 033+4188 and the 
center of Cyg 0B2 obtained by ( lKnodlsedeJl2000bh is 7' 
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The 20-40 keV bandpass considered her e is significantly 
narrower that the 20-60 keV range cited by iDe Becker et alj 
(2007). Based on the power-law spectrum and slope estimate by 
Hsola et al.i (l2007b . we expect the corresponding estimate for the 
20-60 keV range to be a factor 1.8 higher Moreover, the up- 
per limit of the fluxes f or the undetected stars was estimated by 
iDe Becker etall ( l2007h within a PSF of 12'. Using the King pro- 
file describing the project ed spat ial distribution of starlight from 
the association found by (iKnodls eder 2000b), the ratio between 
the flux in the whole association area (radius ~ 60') to that in the 
3EG J20033+4188 region is 6.6. Applying these scaling factors, 
the total luminosity in the whole region of Cygnus should not be 
higher than 7.3 x 10^^ erg s"'. While this scaling is somewhat 
crude, the range of 20-40 keV fluxes we have evaluated is one 
order of magnitude lower than this upper limit. This confirms the 
potential importance of hard X-ray emission from the low-mass 
stars in stellar clusters. 

There are several assumptions that might have led us to 
overestimate the Cygnus OB2 lumino sity. First, the total c lus- 
ter mass of (4-10)x10*Mq derived bv iKnodlsede^ ( 12000 a') has 
been challenged by Hanson (20()3j), who estimated the cluster to 
lie slightly closer at 1 .4 kpc and suggested a total cluster mass 
closer to IO'^Mq. Nev ertheless, based on a study of A-stars in 
the Cygnus OB2 field. iDrewetal] (l2008h find that "a total mass 
of 30 000 - 40 OOOMq would not be surprising". Our adopted 
5 X 10"* Mo is then perhaps up to 1.2-1.7 times too high. 

We have also assumed that all of the observed stellar X-ray 
flux is due to continuous flaring. Some fraction of the observed 
luminosity of these stars might also be attributed to a quiescent 
component. 

One additional factor of uncertainty is related to the dura- 
tion of flares. Indeed, the hard X-rays seen in solar flares that are 
usuaUy associated with the flare impulsive phase generall y decay 
on a more rapid ti mescale than soft X-rays (e.g. Dennis & Zarrd 
ll993 VBenz 2008, and references therein). We considered a ratio 
between the duration of flares in hard and soft band d erived from 
solar measures of flare duration dVeronig et al.l2002l) .The case of 
large stellar flares is not so clear, however. The decay of the 14- 
40 keV flux in the large II Peg flare observed by Swift was not 
obviously shorter than that for the soft X-rays. It should also be 
noted that in large stellar flares there will be a larger thermal con- 
tribution in the 20-40 keV range than in much less energetic so- 
lar flares, in alignment with the general correlation of increasing 
plasm a temperature with flare total energy (e.g. Feldman et al. 
Il995h . Moreover, also the minimum energy of the distribution of 
flares is uncertain, we set it to the median value obtained from 
the analysis of the COUP low mass samples. In our estimation 
we considered that it can be uncertain of two orders of magni- 
tude, that implies a variation of a fac t or ~ 5 in the resulting total 
energy: note, however, that iGuedell (l2009l) found a value even 
lower (few times lO-"*). 

Our assumptions could lead to an overestimation of the real 
value of the total luminosity, but a val ue of ~ 6 ■ \0^^_erg an 
order of magnitude below the De Becker et alj (l2007h flux up- 
per limit, would render the low-mass population of importance 
in searches for emission from specific suspected sources of hard 
X-rays, such as the unidentified y-ray source 3EG 2033-1-4118, 
the unidentified TeV source TeV J2032-I-4130, and massive 
colliding-wind binaries. 

The possibility of detecting hard X-ray emission from nearby 
clustered low-mass pre-main sequence populations provides a 
promising means of investigating high-energy processes on stars 
that are generally too distant to study in detail individually. Of 
particular interest is the relation between soft and hard X-ray 



emission of solar flares (e.g. Ilsola et all 120071 : iBattagha et al.l 
12005 ^. and how the underlying physical basis can be extrapo- 
lated to the other regimes, such as the pre-main sequence stel- 
lar case where hard X-rays from flares can be important agents 
of protoplanetary dis k ionisation (e.g. Ilgea & Glassgoldlll999l : 
iGlassgold et alJl2004 . 

Our discussion assumes that all stars will be "observed" in 
the 20-40 keV range, i.e., that X-rays will penetrate both the 
circumstellar and ambient cluster and line-of-sight extinction. 
Cygnus OB2 is located behind the Great Cygnus Rift that leads 
to extinction of up to (and possibly beyond) Ay ~ 10 (e.g. 
Mas sev & Tho mpson 1991). Extinction of 10"' corresponds to a 
neutral hydrogen column density Nh ~ 1 .5 ■ 10^^ cm^^. The opti- 
cal depth for such a column for 20-40 keV X-rays is completely 
negligible; the Compton ionization and scattering cross-section 
of cTj- ~ 6.65 ■ 10"^^ cm^ renders the interstellar medium opti- 
cally thin to hard X-rays for columns up to Nh ~ iO^* cm^^. 
Except perhaps in rare cases of stars obscured by very dense cir- 
cumstellar disk mid-plane gas and dust, the assumption that we 
will see the whole cluster in hard X-rays should be valid. 

Finally, this consideration raises the possibility of using the 
penetrating power of hard X-rays a s a pointer to our Galaxy's 
hidden superc lusters. iHansonl (l2003h noted that extrapolation of 
the locally-derived Galactic cluster luminosity function indicates 
that our Galaxy hosts "tens to perhaps a hundred" massive clus- 
ters with total mass ~ lO'^M©. These clusters are likely hidden 
behind many magnitudes of extinction and will not be easy to lo- 
cate. Hard X-rays can penetrate such extinction, and the arrival 
of missions in the next decade able to provide relatively precise 
imaging with much greater sensitivity in the ~ 10-100 keV 
bandpass, such as NuSTAR (Harrison et al. 2005) and Symbol-X 
dPareschi & Ferrandol2005b . could provide an assay of this hith- 
erto much overlooked population of our Galaxy. In this context, 
Cygnus OB2 represents a potential "Rosetta Stone", offering a 
nearby super-cluster example that can be well-characterized us- 
ing multi-wavelength techniques that will not be applicable to 
more distant and extinguished clusters. 
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